FALLING OF A VISCOUS JET ONTO A MOVING SURFACE 

A. HLODtt, A.C.T. AARTS 1 ", AA.F. VAN DE VEN+, AND MA. PELETIERt 

Abstract. We analyze the stationary flow of a jet of Newtonian fluid that is drawn by gravity 
onto a moving surface. The situation is modeled by a third-order ODE on a domain of unknown 
length and with an additional integral condition; by solving part of the equation explicitly we can 
reformulate the problem as a first-order ODE, again with an integral constraint. We show that there 
are two flow regimes, and characterize the associated regions in the three-dimensional parameter 
space in terms of an easily calculable quantity. In a qualitative sense the results from the model are 
found to correspond with experimental observations. 
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1. Introduction. In the flow of a viscous fluid jet falling onto a moving surface 
different flow regimes can be distinguished, as is easily observed if one pours syrup 
onto a pancake (Figs 11.11 and ll.2fl . If the syrup is poured from a large height and 
the bottle moves slowly, then the syrup thread is perfectly vertical (Fig. II. 1|) . If the 
bottle is held closer to the pancake and moved relatively fast, however, then the flow 
of syrup becomes curved. 

Note the difference in behavior at the touchdown point. In the straight case, 
the jet hits the surface at right angles, and a little puddle forms on the surface from 
which a thicker 'jet' is transported away by the surface movement (Fig. fTTTf) . Since 
the flow velocity of the surface jet equals the surface velocity, conservation of mass 
implies that at the touchdown point the particle velocity in the free jet is larger than 
the velocity of the surface. In the curved case, however, the jet does not thicken at 
touchdown, and the free jet meets the surface tangentially. 

There is a large body of literature on viscous jets or sheets that impinge upon fixed 
surfaces, where one can observe folding of viscous sheets [?, ?], coiling of viscous jets 
[?] , and viscous fluid buckling [?] . The fact that in these cases the surface is stationary 
is essential; to our knowledge the current work is the first study of a viscous jet that 
falls upon a moving surface. 

In Section we construct a mathematical model of the flow, where we first con- 
centrate on the curved jet of Fie . 11.21 We model the flow as a thin, Newtonian jet of a 
priori unknown length. In Section |3 the original system of equations is transformed 
to a first-order differential equation for the flow velocity and two additional relations 
for two unknown parameters. In Section^] we show that in a certain parameter regime 
the original system admits a unique solution, and we give a convenient characteriza- 
tion of the relevant part of parameter space. In Section we present the solution 
algorithms for the model equations. Results for various model parameters are shown 
in Section [SI and in Section [3 we discuss our results and give some conclusions. 

2. Mathematical model. A thin stream of Newtonian fluid with viscosity r] 
and density p is falling from the nozzle of a bottle onto a moving belt fFig. I2.l|) . We 
use the theory of thin jets (see e.g. [?]) and thus describe the jet as a curve. The 

t Center for Analysis, Scientific computing and Applications, Eindhoven University of Technology, 
Eindhoven, The Netherlands. 

* Author to whom all correspondence should be addressed. E-mail: a.hlod@tue.nl. Postal address: 
Dept. of Mathematics and Computer Science, Technische Universiteit Eindhoven, PO Box 513, 5600 



2 A. HLOD, A.C.T. AARTS, AND AA.F. VAN DE VEN, M.A. PELETIER 




Fig. 1.1. Straight flow of syrup for low surface velocity and large bottle height. 




Fig. 1.2. Curved flow of syrup for high surface velocity and small bottle height. 

magnitude of the flow velocity at the nozzle is v nozz i c , the belt velocity is Vbdt, and 
the distance between the nozzle and the belt is L. The flow of the fluid is stationary 
and the jet has a curved shape. We restrict ourselves to curves under tension and 
therefore require that 

(2.1) Ubclt > ^nozzle- 
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In Lemma 14.11 we show that for a curved jet the flow velocity increases from the 
nozzle to the belt which justifies i|2.1|) . 

The part of the jet between the nozzle and the belt is represented by its center 
line, the curve CD in Fig. 12.21 The point C indicates the nozzle and the point D 
indicates the contact with the belt. The acceleration of gravity is g and the belt 
surface is perpendicular to g. We parameterize the center line by arclength s (s = 
at C, and s — s on( j at D). For each point on the curve we define a local orthonormal 
coordinate system e s (s),e„(s) consisting of the tangent and normal unit vectors at 
the point s. The angle between e s (s) and the belt surface is 0(s). The cross-sectional 
area at s is A(s), and the average velocity of the fluid at this point is v(s) = v(s)e a (s). 




belt 



Fig. 2.2. The geometry of the jet falling from the nozzle onto the moving belt. 



The flow of fluid is described by the equations of conservation of mass and balance 
of momentum [?, (4.18)], which for stationary flow are 



(2.2) 



(A( S )v(s))' = 0, 



(2.3) 



(A{s)v(s)v( S ))' - -(P( S )e 3 (s))' + g A(s), 
P 



respectively, where by ' we denote differentiation with respect to s. The longitudinal 
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force P(s) is obtained from the constitutive law for a Newtonian viscous fluid 

(2.4) P(s) = T)tA(s)v'(s). 

Here r/x is the Trouton elongational viscosity, which for a Newtonian fluid equals 3r] 
[?]• 

Using (|2.2|) and l|2.4|l . we write the balance of momentum (|2.3|) in components in 
the coordinate system e s (s),e n (s), as 

(2.5) As)= 9 -^^ + ,( V -TY 

v(s) \ v[s) 

„x / n^// n gcos0(s) / v'(s)\ , 

where is equal to three times the kinematic viscosity, i.e. fi = 3rj/p. 

Since system H2.5H2.6JI is of second order with respect to the velocity v(s) and of 
first order with respect to the angle 0(s), we need two boundary conditions for v(s) 
and one for O(s). For v(s) we know the velocity of the jet at the nozzle (point C) 
and the velocity at the contact with the belt (point D) 

(2.7) v(0) = w nozzlc , 

(2.8) v(s cnd ) = v hclt . 

Note that the length s en d of the jet CD is unknown. The angle 8(s) at the contact 
with the belt is zero, so 

(2.9) 6(s end ) = 0. 

Because the length of the belt s cn d is unknown in advance, we need an additional 
condition relating s cn d to the distance L between the nozzle and the belt 

/•Send 

(2.10) L= / sin6(s)ds. 

Jo 

The equations l|2.5|l and l|2.6|l together with the three boundary conditions (|2.7H2.9|) . 
and the additional condition l|2.10|l form the complete system for the unknowns u(s), 
6(s) and s on d- 

We next make the equations dimensionless. We scale the length s cn d with respect 
to /tt/ubeit) reverse the direction of s, and move the origin of s to the point D, i.e. 
5 := fbeit(send — s) / fi. The velocity u(s) is scaled with respect to the velocity of the 
belt Ub e it, i.e., Wbeit#(s) := v(s). Also we introduce a new angle 0(s) := 0(s). The 
scaled version of (|2.5H2.10jl reads 

(2.11) U(£) + -^2 =-A- 



(2.12) &(s) = -A 



v(s) J v(s) 
cos0(s) 



(2.13) 5(0) = 1, 

(2.14) W(s en d) = Wnozzle, 

(2.15) 6(0) =0, 

(2.16) / smQ(S)dS = B. 
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Here 

. _ 9H> n _ VbeltL _ fnozzle , _ Send^bclt 

A — — o , i) — , fnozzlc — j and Send — • 

"belt V "bolt H 

All these parameters are positive and B is the Reynolds number. The prime now 
denotes differentiation with respect to s, and in the sequel we omit the tildes. 

3. A first-order differential equation for the velocity. By introducing a 
new variable £(s) 

(3.1) £W=|;W + ^M, 

we can rewrite the equations (|2.11fi2"T2)l as 

(3.2) v(s)C(8) = -AwxQ(s), 

/„ . . — i . . . cos6(sl 

3.3 v(s)e'(s) = -A \> , 

(3.4) v'(s) =£(s)v(s) -v 2 (s). 

For the variable £(s) it is necessary to provide an initial value. To compute it, we 
need to know values of v(s) and v'(s) at the same point. Because we do not know a 
value of v'(s) at any point we prescribe a value of £(s) at s = 0, 

(3.5) f(0) = -Vw, w>0. 

Here we restrict ourselves to a negative initial value for £. Further in this section, see 
<|3.17[1 and l|3.18|l . we explain our choice of the form for the initial value for £(s). The 
value w is unknown in advance and is determined by the requirement that a solution 
of iEOUIPJl has to satisfy the conditions (|!TTo|) and fXTfy . 

Next, we replace the material coordinate s by the time t, according to 

ds = v(t)dt, 

in the system of equations l|3.2H3.4[) together with the conditions (|2.13I^TTB|l and 
H3.5J) . and we obtain 

(3.6) f(t) = -Asine(i), 

(3.7) e'(t) = -aS^M, 

(3.8) t/(t)=£(i) V 2 (i)- V 3 (i), 

(3.9) £(0) = -y^, 

(3.10) 9(0) = 0, 

(3.11) «(0) = 1, 

(3.12) v(t cn d) = Wnozzlc, 

(3.13) [ ^ v(t)sine(t)dt = B. 



Here 



^end 



S ^ d ds 
v(s) 
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represents the dimensionless time necessary to flow from the nozzle to the belt, which 
is unknown in advance. 

To solve the equations i|3.6[l and (|3.7|) we multiply (|3.6[1 by sin 0(f) and (|3.7[1 by 

cos0(f)£(t) and add them, to obtain 

(3.14) (sin0(i)£(i))' = -A. 

We integrate l|3.14f> with respect to t and use the initial condition (|3.10|) to obtain 

(3.15) £(t)8mQ(t) = -At. 

By eliminating sin 0(f) from l|3.15|l and substituting it into (|3.6|> we derive the differ- 
ential equation for £(£) 

(3-16) £'(*) = 



which has the solution 



(3.17) ^) = ±v^+iW- 

Here we have to choose a correct branch of the square root 1(3.17(1 . The branch with 
the positive sign gives negative sin 0(f), see (13.15(1 . which implies an upward-sloping 
jet; the physically reasonable choice is therefore the branch with the negative sign. 
With the initial condition (|3.5|) we get 



(3.18) £(t) = -y/AH*+w, 
and from (|3.15|l and l|3.18|l we find 

At 

(3.19) 0(f) = ar 



^JAH 2 + w 

Summarizing, the problem ((3.6H3.13"|) simplifies to 



(3.20) v'(t) = -v 2 {t){^AH 2 +w + v(t)), 

(3.21) v(0) = 1, 

(3.22) v(tend) = "nozzle, 

(3.23) r Mv[t) — ° 



VAH 2 + 



w 



The unknowns of the problem (|3.2UKl3~2*S|l are the velocity v(t) and the two positive 
parameters w and i C nd- 

4. Existence and uniqueness. We reformulate the problem H3.20hl3~2*3l) as an 

algebraic equation for the parameter w. First we formulate properties of a solution 
v(; w) of for given w > . 

LEMMA 4.1. For any w > 0, equation JtTWJj) has a unique solution v(-,w) : 
[0,oo) -► (0,1] satisfying fQ7)) with v{-,w) € ^([O.oo)). 

In addition, 

1. v(t; w) is a strictly decreasing function oft for fixed w and a strictly decreas- 
ing function of w for fixed t. 
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2. 

(4.1) v(t;w) < 



2 + t^/A T t T 



3. The operator w > v(- ,w) is continuous from [0, oo) to L°°(0,oo). 

Proof. The right-hand side of Ij3.20|l is (7(0) and Lipshitz continuous in v uni- 
formly on fJ, where fl = ({t,v,w} : t G [0, oo),n G (0,1], to > 0). Therefore, locally 
there exists a unique solution of (|3.20|) satisfying (|3.21|) . which continuously depends 
on w [?, Theorem 7.4]. 

From (|3.20(1 it follows that v'(t;w) < whenever v(t;w) > and that v = 
is a solution of this equation. Thus, because of l|3.21f> . v'(t;w) is always negative 
and v(-;w) is strictly decreasing. Since v = is a solution of (|3.20|l . v(t;w) remains 
positive for t > 0. Therefore, v(i; u>) G (0, 1] Vt > 0; this proves the existence and 
uniqueness of v and the monotonicity in t. 

For the monotonicity in w, fix w± > ui2 > 0. Then t/(0;u>i) = — (y / wT+ 1) < 
v'(0,W2) = — (v / ^2 + 1)) an d v(t;u>i) < v(t;u>2) for small t > 0. Suppose that there 
exists a t* > such that v(t*;wi) — v(t*;w2); then v'(t*; w\) > v'(t*; W2), which leads 
to a contradiction with w\ < W2 ■ This completes the proof of part ^ of the Lemma. 

Because v(t; w) > we have 



v'(t;w) < -v(t; w) 2 V AH 2 + w, 



or 



(4.2) f-?T L T > ) > y/AH* + w. 

We integrate (|4.2(l from to t and apply the initial condition v(0; w) = 1 to find the 
following estimate of w): 

, . 2,4 2 
v(t;w) < 7 < ; =. 

2A + WA 2 t 2 + w + w log ( At+V ^^ ) 2 + tV^ 2 ^ 2 + w 

This estimate proves part|21and shows that v(t;w) — > as t — ► 00. 

The right-hand side of l|3.20|l depends continuously on to. This together with the 
estimate <|4.1|) of v(t;w) at i = 00 proves|3 □ 

In order to solve H3.20H3.23[1 we nee d to nn d w f° r which l|3.22|l - 13.23|l are sat- 
isfied. Knowing a correct value of w, we can obtain a solution v(t) which leads to 
a solution of the original problem 1 1|> - (|2 . 16|> . Therefore, next we concentrate on 
finding a correct w. 

Definition 4.2. We define a function I : [0,oo) — > [0, 00) in the following way. 
For given w G [0,oo) let v(-,w) be the solution of \3.2(A\3.21\) given by Lemma \4-l\ 
By items^\ and\^ of Lemma \4-1\ there exists a unique t en d(w) > satisfying 

(4.3) v(t cnd (w);w) = v nozzlc . 
Define I(w) as 

f tc " d(w) Atv(t;w) , 
4.4 I(w) = / ; V ' ; dt. 

Jo y/AH 2 + w 
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By Lemma 14.11 part [3 the integrable function is bounded from above and the integral 
converges. 

Corollary 4.3. Solving \'J.21M3.£3\) is equivalent to finding a w > that satis- 
fies 

(4.5) I(w) = B. 

In the next three lemmas we will show some properties of I(w) which lead to a 
characterization of existence and uniqueness of a solution to (|4.5|) . 

Lemma 4.4. I(w) is a strictly decreasing function ofw. 

Proof. Choose wi and u>2 with 

(4.6) W! > w 2 > 0. 
From part ^ of Lemma 14.11 it follows that 

(4.7) t cnd (wi) < t cnd (w 2 ). 

Combining 14.71) with the statement ^ of Lemma f4. II and (|4.(i|l with the definition of 
I(w), we have 

,W»0^^ ,w<-) Atv{t;Wl) 
v ; Jo vWT^ Jo yJW + m 



rW«>«) Atv(t;w 2 ) 
V 'A 2 t 2 + w 2 



< f 7zr" dt=i( W2 ), 

Jo 



which proves the Lemma. □ 

Lemma 4.5. I(w) is continuous. 
Proof. Fix w > and let 

(4.8) w n — ► w as n — ► oo. 

Then 

/• t - d( "" ) Atv(t;w n ) M 

l(w) — l(Wn) = / =dt— I at 

K 1 V 'Jo \J A 2 t 2 + w Jo VCH 2 + w n 



tendon) 



Jl+J2 



Atv(t; w) Atv(t; w r , 



^JA 2 t 2 + w y/AW + ■ 



■tend(«0 Atv(t;w) , 

dt+ / ; v I ' dt 

'teadKO \/A 2 t 2 +VJ 



Both Ji and J 2 converge to zero as n — * 00; for J\ this follows from the continuity of 
v(t;w) in it) (Xemma I4.H and for J 2 from the continuity of t ea d(w) in w, which we 
prove next. 

From Lemma 14.11 we have that v(-;w) £ C 1 ([0,oo)) and —00 < v t (t;w) < 0. 
Therefore, by the Inverse Function Theorem (e.g. [?, Theorem 9.24]) there exists a 
function t = t(-;w) S (^((0, 1]) such that t(v(t; to)) = t for all t > 0. 

Next note that 

(4.9) V„ := v(tend(w n );Uj) > ^nozzle aS Ti > OO, 

since 

\v{t cnd {w n ); w) - v nozzle \ = \v(t CDd (w n ); w) - v(t cnd (w n ); w n )\ 

< \\v(-;w) - v(-;w n )\\oo — ► 
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by part |3 of Lemma 14.11 Therefore, by continuity of t(-;w) we have 

t C nd(w n ) = t(v(t cnd (w n );w);w) = t(v n ;w) — > t(v nozzlc ; w) = t cnd (w), 

which completes the proof. □ 

Lemma 4.6. lim^^oo I(w) = 0. 

Proof. From the definition of I(w) and v(t; w) £ (0, 1] ( Lemma I4.1|l we have 

i(w) = r M Atv ^ alt < - M alt 



y/w + A 2 t end (w) 2 - ^/w _ At cnd (w) 



2 



A + AH cnd (w) 2 + 

Because t ead (w) decreases in w, by letting w — > oo we find 

lim I(w) = 0. 

i(J— 'CO 

□ 

Summarizing the results of previous lemmas, we formulate a theorem of existence 
and uniqueness of a solution to the original problem (|2.11hl?!TBjl . 

Theorem 4.7. There exists a solution to the problem if^.iitEOT))) if and only if 

(4.10) I(0;A,v nozzic ) > B. 

If it exists, the solution is unique. 

The theorem follows simply from Lemmas 14.41 14.51 and 14.61 




Fig. 4.1. Surface 7(0; A, D nozz i c ). 



In Section we describe two algorithms for computing 1(0; A, w n0 zzic) as a func- 
tion of A and u n0 zzic, resulting in the graph of Fig. 14. II As a consequence of Theorem 
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14.71 a solution to the original problem (|2.11tO?T^|l exists only if the point (A, v nozz \ e , B) 
is below the surface 1(0; A, f nozz i )- 

Note: Non-existence of a curved-jet solution. When the condition l|4.10|l 
is not satisfied, then there exists no solution to the curved-jet equations l|2.11t{0^|) . 
What happens to a viscous jet in this parameter range is not clear. It is possible 
that a stable straight jet exists (9 = it/2); close to the transition between existence 
and non-existence the curved-jet solution becomes straight (see Fig. 16.3(1 . supporting 
this possibility. Other possibilities are coiling (in three dimensions) or buckling (in 
two dimensions), possibly coexisting with a steady straight jet that is dynamically 
unstable. 

At this moment the issue is open, in part since it is unclear which boundary 
condition at the lower end of the jet best reflects the physical situation. We plan to 
return to this question in a future publication. 

5. Numerical approximation. The problem of this paper gives rise to two 
slightly different numerical questions. The first question arises in making a phase 
diagram such as Fig. 16.11 in order to distinguish between existence and non-existence 
of a curved jet we need to calculate 1(0; A, v nozz i e ) and check the existence condition 
1(0; A, w n ozzie) > B (|4.10|l . The second question arises when this condition is fulfilled: 
by Corollary 14.31 we then need to find w > such that I(w; A, i> n0 zzie) — B, from 
which v and can then be determined by solving l(3.20hl3~2"T|l and using l(3.19|l . 

The main differential equation (|3.20|l can be solved either analytically or numer- 
ically, giving rise to two different methods. 

Method 1. When w = 0, it is possible to solve the problem ((3.20ll!OTj) analyti- 
cally (Appendix^). The rescaled domain size z* is then to be determined implicitly 
from 

(2 ,4)i/3 / { J^-J^ VY 1 J_ 8 (V2A/3) 
[ MlOZZlC " Oz*) 2 / 3 [ + [Ji(z*) + J_!(z*) Cl ) ) ' Cl " J f (V2A/3) ' 

where the J a are the Bessel functions of the first kind. We then calculate t e nd(0) and 
1(0; A 

i ^nozzle) 

(6 Z *)V3 J 2 (z*)ci - J_i(z*) 

(^•^) ^cnd(O) = TTR / ( *\ i 7 f~T\ ' 

A 2 ' 3 Ji(z*) + J_i(z*) Cl 



/_ „n w„ . x 1 "^nozzle .^ond(O) 

(5.3) 1(0; A, Wnozzio) = A . 

^nozzle ^ 

This method only is available for the special case w = 0. 

Method 2. Alternatively, one may integrate l|3.20(l numerically until the con- 
dition v(t;0) = Wnozzic is reached. The integral I(w; A,v nozz i c ) can be computed 
numerically as well. This method is available for all w > 0. 

We solve 1(4. 5|l by the bisection method, supplemented with an upper bound on w 
that follows from the estimate ((4.1(1 : since for all w, 
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we have 



^nozzle = v(t cnd (w)]w) < 



^ond(w) 2 



and therefore At en d{w) 2 < 2/v nozz \ e . We thus estimate 

T( , , f uM Atv(t;w) At cnd (w) 2 1 
(5.4) I[w; A, Nozzle) = / , ,„ „ at < — < -=. 

Therefore the solution w of (|4.5|l satisfies the a priori estimate 

1 

W < — o — r. 



nozzle 



6. Results. From Theorem l4.7l it follows that if the parameters A, B, and v nozz \ e 
satisfy l|4.1()|l . then there exists a solution to the stationary curved-jet equations (|2.5I - 
12.1 U|l (or equivalently H3.2UH3.25jl '): otherwise the jet is vertical, or a stationary solution 
does not exist. Fig. 16. II shows the region of existence of such a curved-jet solution. 



s 0.6 




-2 2 4 6 

log A 



Fig. 6.1. Region of existence of a solution to the original problem for B = 1, (grey region). If a 
point (log A, t> nozz i e ) is inside the grey region a curved-jet solution exists (points 1 and 2); in point 
3 no such solution exists, and we conjecture that an actual jet is straight. If a point (log A, « nozz i e ) 
is closer to the border of the grey region a shape of the jet is more vertical (point 2). 



As a reference configuration for the numerical experiments shown below we con- 
sider syrup with viscosity 77 = 3.2 Pas and density p — 1000 kg/m 3 pouring from 
the height L = 2 ■ 10~ 2 m. The velocities of the belt and the flow at the nozzle are 
Vbelt = 0.5 m/s and u n ozzic = 0.05 m/s, respectively. 

Fig. 16. 2l shows curves in non-dimensional parameter space (A, v nozz i c , B) corre- 
sponding to variation of a single (dimensional) physical parameter L, fx, fbcit, or 
^nozzle- In the figure, we see that if fj, or Ubdt decreases, or if L or v llozz \ c increases, the 
point (A, Wnozzic, B) eventually leaves the region {(A, v nozz i c , B) : 1(0; A, u noZ zie) > B}. 
Close to this transition the curved jet becomes vertical. In Fig. 16.31 we present the 
shapes of the jet for specific values of the parameters along each of these curves. 

Summarizing the numerical experiments, we observe that by increasing the flow 
velocity at the nozzle or the distance between the belt and the nozzle the jet shape 
becomes more vertical; the same is true if we decrease the velocity of the belt or 
the kinematic viscosity. The jet becomes exactly vertical when the parameter point 
(A, fnozzie, B) approaches the critical surface {I (0; A, v nozz i e ) = B}. 
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B 



2 




6 



0.4 V, 



0.8 



nozzle 



log A 



Fig. 6.2. Curves in non-dimensional parameter space (A, f nozz i e , B) as we change one of the 
process parameters (L, fi, «belt , v nozzle)- The grey parts of the curves below the surface 1(0; A, f n0 zzle) 
correspond to the curved jet; we conjecture that the black parts of the traces correspond to a vertical 
jet. Line a: increasing L; line b: decreasing pi; line c: decreasing i>b e it>' and line d: increasing v nozz i e . 

7. Conclusions. In this paper we propose a mathematical model of the falling 
of a viscous jet onto a moving surface. We assume that the jet is falling under gravity 
and has a curved shape. The model consists of two differential equations, one for the 
flow velocity and one for the angle describing the jet's shape. An additional relation 
fixes the unknown length of the jet. 

The initial system of equations is partially solved and then transformed to a 
first-order differential equation for the velocity. By introducing an additional scalar 
parameter w the problem is reformulated as an algebraic equation for w (|4.5() . For this 
equation we formulate an existence condition (|4.10() and prove uniqueness, thus giving 
a complete characterization of existence and uniqueness for the original equations. 
Finally, we solve the equation for w numerically and recover the solution of the original 
problem. 

We have shown that if the existence condition (|4.10(1 is satisfied, then the shape 
of the jet is curved; we conjecture that in the alternative case the jet is vertical, 
but this case lies outside of the scope of this paper. Furthermore, the model shows 
that the curved jet becomes more vertical when: i) the distance between the nozzle 
and the surface increases, ii) the flow velocity at the nozzle increases, iii) the surface 
velocity decreases, or iv) the kinematic viscosity of the fluid decreases. These results 
correspond with those observed in the basic experiment described in the introduction. 

Appendix. Calculation of 1(0; A,v nozz \ e ). First we calculate v(t;0) analyti- 
cally. The differential equation for v(t;0) follows from (|3.20() and l|3.21|l . 



(A.l) 



u'(t;0) 



v 2 (t; 0)(At + v(t; 0)), v(0; 0) = 1. 



By replacing v(t; 0) by Z(t) = l/v(t; 0), we find 



(A.2) 



Z'(t)Z(t) = AZ(t)t + 1, Z(0) = 1. 
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y(cm) 



Line 


L(cm) 


1 


2.00 


2 


2.46 


3 


2.74 


4 


2.84 




0.1 0.3 0.5 0.7 0.9 1.1 1.3 



x (cm) 




Line 


H (m 2 /s) 


1 


0.0096 


2 


0.0067 


3 


0.0059 


4 


0.0054 



0.1 0.3 0.5 0.7 0.9 1.1 1.3 



x (cm) 



(a) Varying L. 



(b) Varying fi. 




Line 


V b .i,(m/s) 


1 


0.50 


2 


0.28 


3 


0.25 


4 


0.23 



0.1 0.3 0.5 0.7 0.9 1.1 1.3 




Line 


Vnozzle (m/s) 


1 


0.05 


2 


0.10 


3 


0.12 


4 


0.14 



0.1 0.3 0.5 0.7 0.9 1.1 1.3 



(c) Varying ?j bclt . 



(d) Varying u noZ zle- 



Fig. 6.3. The shapes of a jet for different values of the process parameters (L, /j, Ubelti ^nozzle)- 
The reference values for the parameters are rj = 3.2Pas, p = 1000 kg/m 3 , L = 2 ■ 10 -2 m, t>bclt = 
0.5 m/s, and f n0 zzlo = 0.05 m/s. 



We seek for a solution of 1A.2|) in parametric form. With the substitution 

(A.3) Z{z) = z + A/2t 2 {z), 

where z is a parameter, (|A.2J| becomes 

(A.4) = At 2 (z) + z, t(l)=0. 

Here the initial condition is deduced from (|A.3J| by setting = and Z(z) = 1. 
This differential equation is known as the special Riccati equation [?, p. 4, type 4] 
and has the solution 

(A.5) t(z) - V a V 3 y - 



with 



J_i(a) , 

(A.6) Cl = ^^ ; a=y/2A/3. 

Here the functions J a are the Bessel functions of the first kind. The velocity v is given 
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by 



(A.7) v{z ' 



To write the result in a more elegant form we replace the parameter z by z 
V2Az 3 / 2 /3 (we then omit tildes) 

,, D , , , (6Z) 1 / 3 Mz)ci~ ■/-§(«) 

(A ' 8) ^ = 4(z) + J_,(z) Cl - 

(A - 9) B(z)= kF[ 1+ ( fefe y 

To calculate i e nd(0) from the solution (|A.8|I and (|A.9|I it is necessary to find z* satis- 
fying 

(A ' 10) ^-(3^)2/3 1+ U(z*) + J ,(z*) Cl 



and then substitute z = z* into (|A.8|1 . 

The equation (|A.10|I has many solutions. A correct solution z* is the first solution 
of ijA.lOfl after the point a. It is convenient to search for z* in the interval (a, zq) 
using the bisection method [?]. Here, z$ is the first zero of v(z) according to (|A.9(1 
after the point a. 

Next, we have to find a correct zq. Because zeros of v(z) coincide with zeros of 
(A.ll) Ji(» + J_i(z)ci, 

we can look for the first zero of IjA.llf) after a. Using l|A.6|) . we can rewrite the latter 
as 

(A.12) J i (z o )J i (a)+J_ i (z )J_ i (a)=0. 

This equation can be rewritten in terms of Airy functions [?, 10.4.22 and 10.4.27] as 

2 2 

(A.13) Bi(-z )Aif(-a) - Ai(-z )Bi'(-a) =0,z = (^j ' , a = (^j ' . 

Using the representation of Airy functions via modulus and phase [?, 10.4.69 and 
10.4.70] 

Ai(-z ) = M(z ) cos 6 (z ), Bi(-z ) = M(z )sm9(z Q ), 



Ai'(-a) = N(a) cos 0(d), Bi'{-a) = N(a) sin 0(d), 
we see that (|A.13|I becomes 

(A.14) sin(0(z o ) - 0(d)) = 0. 
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For large z ^> 1 and a ^> 1 the asymptotic expressions for 9(z ) and (j)(a) [?, 10.4.79 
and 10.4.81] are given by 

M» \ 71 3 ~ 2 / 3 fi 5 , 1105 



and 



4 2 u V 32 ^o 6144zg 



, N 3tt 3 ,,o / 7 1463 



4 2 V 32d 3 6144a 6 
or in terms of zq and a <|A.13|) 

1105 



(A - 15) 4 "* I 1 "tS? "31104* 

and 

/4 , . N 3tt / 7 1463 \ 

^^T-H 1 + ^"3U04^J- 

After substituting (|A.15|I and (|A.15f) into l|A.14l) for a > 1 we find 
(A.17) z ~a + 7r/2. 

When a is not large we can find zq numerically by looking for a solution of l|A.12|) in 
the interval (a, a + n). 

Once zq is found we find z* and consequently compute i cn d(0). Knowing £ cn d(0) 
we can compute 1(0). To avoid computation of the integral 

/•tend(Q) 

7(0) = / v(t;0)dt, 



we can calculate this integral using the differential equation (jA.lfl . when written as 



1 w 



(A.18) =^ + «(*;0). 

By integrating this equation from to i O nd(0), we get 

(A.19) ( 1 —— - —^—) = A tcnd( ° )2 + / tC " d(0) v(t; 0) dt. 

V ' v«(W(0);0) v(Q;0)J 2 7 1 ' ' 

We use the definitions of 1(0) and i e nd(0), together with the initial condition w(0; 0) 
1 to obtain 

(A.20) 1(0) = j^li^i _ A !ggd(0)^ 

^nozzle 2 



